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A recent experiment reports a non-local spin-signal that shows oscillatory behavior as a function 
of gate voltage when the contacts are magnetized along the direction of current flow, but not 
when they are magnetized perpendicular to the current, in agreement with the predictions from a 
simple theory. In this paper we first present a straightforward extension of this theory to include 
the angular spectrum of electrons and the extended injecting and detecting contacts. The results 
are in good agreement with those from a non-equilibrium Green function (NEGF)- based model 
with contact parameters adjusted to fit the experimental contact conductances. They also describe 
certain aspects of the experiment fairly well, but other aspects deserve further investigation. 

PACS numbers: 



I. INTRODUCTION 

Voltage controlled spin precessionfsee [l[ and refer- 
ences therein), proposed in 1990 posed two diffi- 
cult challenges namely (1) spin-polarized injection into 
a semiconducting channel, and (2) gate control of the 
Rashba spin-orbit interaction (RSO) in the channel 
The latter was demonstrated by Nitta et al. in 1997 us- 
ing an inverted InGaAs/InAlAs quantum well with a top 
gate Q. But spin-polarized injection into a semiconduc- 
tor proved to be a more difficult challenge which has only 
recently been overcome through the combined efforts of 
many groups around the world Q. Very recently, Koo et 
al. [l| have combined both ingredients, spin-polarized in- 
jection and gate-controlled RSO, into a single experimen- 
tal structure using a high mobility InAs heterostructure 
with a top gate interposed between the current contacts 
and the voltage contacts (Fig. [TJa)). The non-local volt- 
age signal showed an oscillatory behavior when the 
contacts are magnetized along the direction of current 
flow, but not when they are magnetized perpendicular to 
the current flow (Fig. lUb)), as expected from the theory 
presented in [2[. Furthermore, it was shown in [l[ that 
the oscillation (see Fig. [T](b)) is described well by the 
expression 

K.p = Acos( + ) (1) 

where m* is the effective mass, a{VG) is the RSO mea- 
sured independently from the Shubnikov-de Haas (SDH) 
beating pattern, and A and 4> are fitting parameters. The 
oscillation period 2m*a{VG)L/h^ was derived by Datta 
and Das Q for electrons with wavevectors that are purely 
along the direction of current flow (x-) with fcy = 0, not- 
ing that 'in practice we have an angular spectrum of elec- 
trons' and the 'effect is reduced as k turns away from the 
X-axis.' 
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In this paper we will first (Section [IT]) describe a 
straightforward extension of the theory in to in- 
clude the sum over the angular spectrum ky of electrons. 
The results closely follow those obtained earlier in ref. Q 
which are more general since they include both electron 
and hole systems. We show that the results from this sim- 
ple model are in good agreement with those in (FiglTJd)) 
from a non-equilibrium Green function (NEGF) based 
model with contact parameters adjusted to fit the ex- 
perimental contact conductances fSection lIII[) . We hope 
that a careful comparison of experiments will help re- 
fine model proposed here and establish this effect on a 
firm footing, so that it can be used both for fundamen- 
tal studies as well as for various proposed applications 
such as spin-filtering, magnetic recording and sensing or 
quantum computing Q. 

II. SIMPLE MODEL 

We start from an effective mass Hamiltonian for a two- 
dimensional conductor with RSO interaction of the form 
(a: Pauli spin matrices) 

H - + ^^''^^y - '^^^^) (2) 

We neglect Dresselhaus spin-orbit (DSO) coupling since 
this is believed to be small in structures of this type 0] . 
Eq. [2] leads to the dispersion relation 

E = ^ ±ak,k = +^kl + kl (3) 

with the upper and lower signs corresponding to eigcn- 
spinors of the form {1 ± exp(i(/))}"^, where tan0 = 
— kx/ky- Here, x and y are the longitudinal (or trans- 
port) and transverse direction respectively following the 
co-ordinate system used in [l[, which is different from 
that used in [2| . Assuming periodic boundary conditions 
in the transverse direction, both E and ky are conserved 
in the absence of scattering and the two eigenmodes have 
different kx's so as to satisfy Eq. [3] with the upper and 
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FIG. 1: (a) Schematic structure and (b) experimentUy ob- 
served non-local signal in [11]. Calculated non-local signal 
for the structure in (a) using: (c) simple analytical model 
(section |ll]) and (d) using NEGF based model (section . 



Parameters: Pc = 6.8%= 
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^"T^" and Gm = 4 X 

W'"/7n^/Q (Gm( 

being the contact conductance per unit 
area for majority (minority) spins), carrier density ns = 2.7 x 
lO^^m"^ FM contact lengths La = 0.2^tm, Lcd = 0.25^im, 
FM contact spacing Lch ~ 1.65/im, width W — 8/im, effec- 
tive mass m* — 0.05mo. 



lower signs respectively. For small a we can write ap- 
proximately 



2m*a 



(4) 



with fco = V2m*E/h. 



We define the spin- voltage Vx (or V^) as the differ- 
ence between the voltages measured in the parallel and 
the anti- parallel configurations by X- (or Y-) directed in- 
jecting and detecting magnets. This is expected to be 
twice that measured in the parallel configuration using 
a setup like the one shown in Fig. [TJa) flO| which is ex- 
actly how the oscillatory signals Vx in [l[ are measured, 
although the spin- valve signal Vy is measured the same 
way our Vy is defined. It is shown in Appendix A (sup- 
plementary information) that for a point injector located 
at a; = and a point detector located at a; = L, as- 
suming ballistic transport, the voltage signals for X- and 
Y-directed magnets can be written as {B: constant) 



(5a) 



Vxo{E,ky) = B + (1 - s^) cos 
Vyo{E, ky) = B s^) + cos 



(5b) 



where s = ky/ko = hky /\/2m*E and 9l = 2m*aL/h?. 
The contributions from different ky all act 'in parallel' 
giving a voltage equal to the average. At low tempera- 
tures we can average the contributions from all transverse 
wave- vectors ky over the Fermi circle {E = Ep) to write 



X(Y) 



(Mvy 

-fco 27rfco 



{EF,ky 



(6) 



Interestingly, the results obtained from the integration 
in Eq. [5] looks almost exactly like the single cosine re- 
sult in Eq. [T] that describes the experimental observa- 
tions. This can be understood by noting that the argu- 
ment ^i/Vl — has a stationary point at s = [11| 
and we can use the method of stationary phase to write 
approximately. 
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As shown in Appendix B (supplementary information) 
these approximations describe the results from the exact 
integration quite well for 0^ > 2-k which is true for the 
range of a and L involved in the experiment. This should 
help answer some of the concerns raised in a recent com- 
ment [ill. 

Let us note that the simple results presented above are 
made possible by our assumption of periodic boundary 
conditions (PBC) in the y-direction making ky a 'good 
quantum number' like E. Most of the prior work, on 
this topic [i3|, uses hardwall boundary condition (HBC) 
which does not seem to permit the simple decoupling of 
different transverse wavevectors (fcy) due to non-trivial 
'boundary scattering'. We have checked numerically that 
the use of HBC does not change the conclusions described 
above in a significant way although some details are dif- 
ferent. Furthermore, one could argue that since the ac- 
tual boundaries in the experimental structure of [l[ are 
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relatively far away (W = Sfim) the physics is better cap- 
tured by a model employing PBC like ours. However, 
the possible role of boundary scattering deserves further 
attention. 

It is interesting to note the similarities and differences 
between this simple model for the voltage controlled spin- 
precession signal and the usual model for the Hanle sig- 
nal H 



Vh 



at — , cos — - — I e ' ^ 



(8) 



The cosine functions in Eqs. [5] can be written as 
cos(2akQt / h) where t is the transit time L/vx = 
m*L/hkQ\/\ — s^, showing that 2afco in our problem 
plays the role that gfisB plays in the Hanle precession 
signal. Since we assume ballistic rather than diffusive 
transport we have a different weighting function for dif- 
ferent transit times. But the most important difference 
is that Hanle signals are typically observed within tens of 
Gauss around B = while the experiment we are analyz- 
ing has a varying between 8 x 10~^^ and 12 x 10~^^eV- 
m with /co = 4.1 X 10*m~^ corresponding to values of 
\g\B close to ^ 140T as noted in [l|. At such high val- 
ues of \g\B, the Hanle signal is usually reduced essen- 
tially to zero because of the spread in the transit time 

caused by diffusive transport. One would expect the 
same in the present case, were it not for ballistic trans- 
port. By contrast, the Hanle signal around B = is 
relatively robust and it would be interesting to look for 
an analogous voltage-controlled signal in shorter struc- 
tures or perhaps in structures where the RSO, q;(Vg), 
can be tuned through a = [l^ . 

The simple model here makes no prediction about the 
amplitude B, but it does suggest that the peak-to-peak 
amplitude of the oscillation in Vx should be 3tt/^/2tt6l 
times the spin- valve signal Vy- Using m* = O.OStoq: ~ 
10~^^eV-m, L = 1.65/im, this suggests Vy = 1.2Vx(p- 
p). Experimentally, the p-p oscillatory signal ~ GfiY 
which equals Vjs:(p-p)/2(since the experiment measures 
the parallel-antiparallel difference we have defined as Vx) 
but Vy is only ~ 6/xV. Possible reasons for the discrep- 
ancy are discussed at the end of this paper, but here we 
would like to note that we expect a further reduction in 
the amplitude of the oscillatory component due to the ex- 
tension of the injecting and the detecting contacts along 
X giving rise to a spread in the values of 0l in Eqs. I5al 
and lSbl We can write 

Bcos{eo + 9, + 9d + Tr/A) 



Vx = CiCd- 



(9) 



where Ci and Cd are numbers less than one represent- 
ing the averaging effects of the injecting and detecting 
contacts respectively and 9i , 9d are the additional phase- 
shifts introduced by the injecting and detecting contacts 
respectively in addition to ^o- To estimate 9i, 9d or 
Ci, Cd we need to know (1) the spatial uniformity of 
the injecting and detecting contacts, (2) how the elec- 
tronic wavefunction evolves under the contacts, and (3) 



how the RSO a{VG) varies under the contacts. Re- 
garding point 1 we assume the contacts to be uniform 
and the NEGF model described next should account for 
point 2 within this assumption. However, point 3 re- 
quires a careful treatment beyond the scope of this pa- 
per. Here we simply note that Eq. |9] describes the shape 
of the oscillatory Vx{Vg) quite well with the following 
choice: = 'i'^n* a{yG)Lc h / h? , where Lch = 1.65^m 
which is the experimental center-to-center distance be- 
tween contacts and 9i^d = m*a{VG = 0)Lci^d/fi^ where 
Lci = 0.2^m and Led = 0.25^m equal to half the con- 
tact widths. The result from Eq.[5]also matches that from 
the NEGF model (see Fig. [Ijd)) to be described next in 
shape and amplitude if we use d^d = sm{9i^d) / Oi^d which 
can be justified if the electronic wavefunction is assumed 
to remain constant under each contact. 



III. QUANTITATIVE NEGF BASED MODEL 

One way to make the results from the quantum 
transport model quantitative is to use the NEGF-based 
method described in detail in [l6j . The inputs to this 
model are the Hamiltonian [H] and the self-energy ma- 
trices [E] (Fig [2]). For H we use a discrete version of 
the one used in Section [H] (Eq. [2]), as described in [l^ 
assuming PBC along y as discussed above. We ne- 
glect all scattering processes since both the mean free 
path and the spin-coherence length are believed to be 
longer than the longitudinal dimensions at low temper- 
atures. To understand the signal decay at higher tem- 
peratures will require a consideration of both momen- 
tum and spin-relaxation processes, but we leave this 
for future work. The self-energies for the FM contacts 
have the form — (i/2)7[/ -I- Pca.h] where the 
polarization, Pc = {Cm — Cm) /{Cm + Cm) and n is 
the unit vector in the direction of the magnet. The 
constant 7 = 'k{Cm + Gm)?i'^/e^m* is chosen to give a 
tunneling conductance equal to the experimental value 
. The NM contacts (Ei,E4) are represented similarly 
with Pc ~ 0. Finally, the long extended regions outside 
the channel at two ends (see Fig. [Ija)) are represented 
by two semi-infinite contacts whose coupling is given by, 
^^(7?) = '''L{R)gsT^j^^ where t is the spin-dependent cou- 
pling matrix between the contact and the channel and gs 
is the surface Green's function. The transmission func- 
tions are calculated from the NEGF model and contacts 
3,4,L and R are treated as voltage probes with zero cur- 
rent (following the approach introduced by Buttiker, see 
section 9.4, in [13 )■ Note that although we are not in- 
cluding scattering processes explicitly, the voltage probes 
introduce an effective spin-scattering that reduces the 
signal. Indeed both Vx and Vy increase significantly if 
we remove the end regions represented by and E^;. 

For contacts 1 and 2 we adjust the applied potential 
difference (fii — /i2) to obtain a current level equal to the 
experimental value. The voltage signal ii^p is obtained 
from the difference between the /i3p measured with par- 
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FIG. 2: NEGF based lateral transport model for the struc- 
ture in Fig. [TJa), with E2 and E3 representing injecting and 
detecting ferromagnetic contacts, Ei and E4 representing as 
non-magnetic contacts (NM), E_t and E/j representing the 
semi-infinite regions outside the central region 



allel contacts and fisAP measured with anti-parallel con- 
tacts. We use a contact conductance of Gm + Gm = 
4 X lO^^fi^^m"^ based on the experimental parameters 
in Q and a Pc = (Gm - G^)/{Gm + G,„) = 0.068 to 
match the spin- valve signal, Vy- Fig. [Ijd) shows the nu- 
merical results obtained using a channel with a(VG') of 
length LcH = 1.65/^m and contacts with fixed a{VG — 0) 
of lengths Lqi = 0.2/xm and Lcd = 0.25/zm. The os- 
cillatory signal matches the experimental observation in 
shape but the amplitude is smaller. One possibility is 
that the Pq we use has been calibrated for the spin- 
valve signals obtained with Y-directcd magnets. The 
same magnets when forced into the X-direction for the 
oscillatory signals may have a higher effective Pc espe- 
cially since no anti-parallel measurements are involved. 
However, to account fully for the discrepancy a signifi- 
cant increase in Pc would be needed and other sources 
of discrepancy should be investigated. 



IV. DISCUSSION 

In summary, we have presented (1) a straightforward 
extension of the Datta-Das theory Q to include the 
angular spectrum of electrons and the extended con- 
tacts, and (2) a more elaborate model that treats the 
actual non-local experimental structure using a NEGF 
based approach widely used in nanoelectronics. The sim- 
pler theory provides a number of insights and is well- 
supported by the more elaborate model, identifying sev- 
eral features that deserve further investigation. Specifi- 
cally it seems that while the experimental oscillation pe- 
riod shows good agreement with theory, the amplitude 
relative to the spin-valve signal Vy is larger showing es- 
sentially no reduction expected from averaging over the 
angular spectrum and the extended contacts. Possible 
reasons for this discrepancy deserve carefull attention. 
On the theoretical side, it is possible that the contribu- 
tion from high ky components are suppressed because 
they have shorter effective spin coherence lengths and a 
purely ballistic theory misses this aspect. Another pos- 
sibility is that the effective polarization Pc is lower in 
the antiparallel case which affects only Vy and not Vx- 
The structure and nature of the injecting and detecting 
contacts also require careful consederation. 
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Supplementary information 



Appendix A: Derivation of Eqs. 5a and 5b 



We start from (see Eq. 3) 



E 



2m" 
2m' 



■{k\_ + kl.) ~ a J k\^ + fcf. 



and 



1 

+ cxp(i(? 



1 

cxp(i^ 



(Al) 



(A2) 



where tan0_|_(_) = — ' . The incident state {?/';} is written a hnear combination of {^/'+} and {'0-} 



(A3) 



where [4*] = [{?/'+}{?/;_}]. After propagating from x = Q to x = L, the final state is written as (^+(-) = kx+(-)L) 



{iPf} = Acxp(i6l+){^A+} + Bcxp(i6l_){V'-} 



= m 

Hence we can write, {ip/} = with 

M = m 

where, 



exp(i0+) 
exp(i6'-) 



cxp(i0+) 
exp(i6'_) 



1 1 

cxp(i0+) — exp(i(/)_ 



(A4) 



(A5) 



(A6) 



Multiplying out the matrices leads to 

Me 



exp(i(/)+ + i6_) + exp(i0_ + i6^) exp(i6'+) — exp(i0_) 

{cxp{i0^) — exp(i6'_)} exp(i0+ + i(f)-) cxp(i0+ + + exp(i(/)_ + 



(A7) 



exp(i(/)+) + cxp(i(/)_) 

Setting ^ (j>- = (j> (this amounts to ignoring non-orthogonality of "+" and "-" states), the expression simplifies to, 



[t] 



Note that [t] can also be written as 



exp(z0+) + cxp(i9^) {exp(i0_|-) — cxp(i0_)} exp(— iij 

{exp(i6'+) — exp(i0_)} exp(i0) exp(i6'_|-) +exp(z6'_) 
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(A8) 



(A9) 



where A9 = 9^ — 9_ = — and fi is a unit vector in the direction of the effective magnetic field: n = 

cos0i + sin0y. This form is intuitively appealing showing the transmission [t] as a product of a simple phase-shift 
exp {i (9^ + 9 ) /2} and a rotation around n by A^. 



For z-polarized contacts in parallel configuration, 
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and in anti-parallel configuration. 
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For y-polarizcd contacts in parallel configuration, 
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and in anti-parallel configuration. 
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Noting that taut/) « —kx/ky and fcg « fc^ + ky we can write. 
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Appendix B: Derivation of Eqs. 7a and 7b 




FIG. 3: Left: Numerical calculation (squares) of Eqs. 6, 5a versus analytical expression (solid) in Eq. 7a as a function of 
Right: Numerical calculation (squares) of Eqs. 6, 5b versus analytical expression (solid) in Eq. 7b as a function of a. 



From Eqs. 6 and 5a, 



Vx = 
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Noting that the phase has a stationary point at s = we expand it in Taylor's series around s = to obtain 
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as stated in Eq. 7a. 

Similarly from Eqs. 6 and 5b, 
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B B TT, 
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s^ exp 



as stated in Eq. 7b. 



